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PREFACE 



Precipitation static has long been recognized as a particularly serious 

problem in aviation radio navigation and communication systems. Many investi- 

8 , 11 , 20 

gations as to its causes and nature have been made with successful 

progress, especially during the past fifteen years. Static free navigational 
systems (VOR Omni-range), for short ranges, and line of sight communications 
systems (VHF-UHF), are now standard equipment on all aircraft, although it 
must be pointed out that these solutions involve a sort of strategic re- 
treat, that is, a shift to higher frequency levels where the precipitation 
static energy content is of a lesser value. Also the investigations that 
were made all utilized propeller driven aircraft for their test vehicles 
and although projects were earmarked to include jet type of aircraft, none, 
to the author's knowledge, were ever reported on. 

The problem that is the subject of this thesis was essentially to de- 
sign a device to measure the relative noise power or voltage due to corona 
discharge, at certain specified locations in the frequency spectrum, and 
be carried aboard jet type aircraft* This is a bit of a swift explanation 
but a detailed account will be given later. 

Also an attempt is made to correlate information regarding some of 
the statistical aspects of noise analysis and although it is not a completely 
rigid treatment, it describes in a condensed manner, a description of noise 
signals and their effects on electronics circuits. 

The author would like to express great appreciation to Messrs. H. P. 
Blanchard and Maurice Mills of the Stanford Research Institute for 
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their infinite patience and guidance, and to Professors Robert Bauer and 
Robert Kahal for encouragement and interest# 
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Chapter I 
INTRODUCTION 



1. History and General Considerations of the Problem 

Precipitation static may be described as radio interference caused by 
corona discharge. It is sometimes labled as, 'Poisson noise', is tech- 
nically named, 'impulse noise', and is statistically described as, 'nearly 
normal random noise'. It is but a single type of the many kinds of noise 
that have been studied, analyzed and catalogued, although it does have the 
dubious distinction of being the number one culprit of aviation naviga- 
tional and communications systems. 

The first analysis and explanation of precipitation static was that 
charged particles, created by the turbulence and charged energy of thunder 
clouds, generated static when striking the antenna. No attempts were made 
to alleviate these effects for many years until the introduction of the 
shielded loop antenna in 1931*. It was assumed that since particles did not 
strike the antenna proper, static would be reduced. However, this assump- 
tion was not substantiated in practice because it was found that in actual 
air operation, severe static easily disabled the system. 

In 1939 after more investigation, a new theory H assigned the main 
cause of precipitation static in aircraft to that of corona discharge. 

Later the concentrated effort of the Army-Navy Precipitation Static Project ® 
verified this theory which resulted in the subsequent changes to air-frame 
construction and airborne electronic equipment that are in standard prac- 
tice today. 

With the increased use of jet and rocket propulsion in supersonic 
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aircraft and guided missiles, it was erroneously assumed that the determi- 
nental effects due to precipitation static would be of no consequence, 
since in this type of craft, a natural discharging process due to the ther- 
mal ionization of the exhaust gases, should efficiently transfer free charge 
to the atmosphere. This is partially accomplished because a complete recom- 
bination of the ions produced in burning the fuel has not as yet taken place, 
thus the electric field of the aircraft will capture those ions of opposite 

O 

sign while simultaneously repelling those of the same sign . 0 

In a recent report by R. L. Tanner, ^it was shown that the strength 
of the noise signals due to corona discharge is directly related to the 
dipole moment of the discharge current . The dipole moment is proportional 
to the product of the current flow and the distance over which it takes 
place. Thus in any supersonic aircraft which strikes ice crystals (cumulo- 
cirrus type clouds), or heavy turbulence (cumulo-nimbus type clouds), or 
the layers of ionized air located at extremely high altitudes, a static charge 
builds up in the vicinity of plastic antenna housings, canopies and control 
surfaces. It is here that one is faced with a new type phenomenon that is 
based on a very old principle: that substances of high dielectric constant 
acquire a positive charge when rubbed with substances of lower dielectric 
constant. This is called triboelectric charging, and can be visualized as 
a scrapping off of the electrons by the airstream. Naturally this charging 
eventually culminates in a corona type discharge at many points close to 
the outer perimeter of plastic housings and the dipole moment of the cur- 
rent discharge acts as an efficient antenna to couple this electrical 
noise into the receiver input circuits. 

The problem specifically, is to determine the average noise power or 
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the average voltage produced by this static, recorded for a period of 
time equal to the flight duration of the aircraft. This system must dupli- 
cate an actual navigational system's characteristics with respect to band- 
width, resonant frequency, servo system time constants and antenna match- 
ing stage. 

The actual navigational system can be represented in general by the 
diagram in Figure 1. 
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TIME CONSTANT 
RC = 100 seconds 



FIGURE 1 

Simplified Diagram of Navigational System 
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To illustrate: 



If the incoming waveform is sinusoidal in nature, then the output of 
the mixture will be the sum and difference terms of the original inputs. 
This is applied to the integrator which acts to average only the difference 
terms since it is a low pass filter. The output will then be a very slowly 
varying voltage and the subsequent error signals will also be very small 
in magnitude. 

Where precipitation static is present the waveforms (idealized) might 
take the form of those shown in Figure 2. Since power is proportional to 
the voltage applied to the servo system, the disturbance to the system is 
related to all the input noise pulses integrated over a period of time 
equal to the time constant of the system, and can be written as 



The error signal to the reinsertion oscillator will then be proportional 
to the input average noise power so that periods of heavy static would pre- 
cede large error signals. 

The system that was designed and constructed duplicating the above 
characteristics is now being operationally used to record the average noise 
power over periods of ten to twelve hours. These results are later com- 
pared to the navigational errors plotted separately by the navigator to see 
what relative values of average noise power cause the maximum errors in 
the system. 




Phase Disturbance in System 
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(A) 












FIGURE 2 

Idealized Waveforms 
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(D) 








FIGURE 2 

Idealized Waveforms 
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A second system to measure the average noise power at specific fre 
quencies is also described in Chapter III with some data which shows re 
suits of actual flights made using the system. 
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2. The Projected Method of Solution: 

A study of the available systems for measuring noise due to atmospheric 
sources and corona discharge indicates a common tendency to make the major- 
ity of measurements using an electric field intensity meter, or circuits of 
some such similar nature. The results are then interpreted using a graph- 
ical portrayal of corona discharge current as a function of electrical field 
intensity. In this problem it is desired to measure the energy components 
of the noise at a specified frequency and known bandwidth, therefore a 
relationship with respect to time must be set up to measure this noise power. 
But the instantaneous noise voltage induced in the antenna cannot be con- 
sidered as a regular function of time, for if it were, both amplitude and 
phase relationships would have to be known for one complete period. Since 
it cannot, the problem is one that is statistical in nature. 

Statistics may be succinctly defined as the science of reduction and 
analysis of data. Its application to the problems associated with com- 
munication processes can be mainly attributed to N . Wiener^ who was one 
of the foremost contributors during the development years of statistical 
prediction theory, particularly concerning stationary random time func- 
tions, a random time function being defined as a function which does not 
follow any specific law of variation and therefore is not precisely predic- 
table. The type of noise being investigated in this paper falls into this 
category for it consists of a superposition of impulses occurring at ran- 
dom times with overlapping individual pulses. It is sometime referred to 
as Poisson noise because the Poisson distribution (Appendix I) is used to 
describe its statistical structure. In statistical applications the 
Poisson Distribution often appears when concerned with the number of 
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occurrances of a certain event for a very large number of observations. 



3,5 



where probability for the event to occur in each observation is very small. 

It is for this reason that it is sometimes called the "small probability law". 

For the case of very hij^i density Poisson noise, that is, when the number 

of elementary random impulses is large, and their mean duration is not 

too brief, a heavy overlapping occurs and the result is nearly normal ran- 

17 

dora noise. This, according to Middleton 1 can be considered to have a 
normal distribution with negligible correction terms, and the shape of the 
pulse is unimportant. (Appendix I) 

It is plausible to assume that the corona discharge that takes place 
will have the characteristics of high density Poisson noise, for when the 
same navigational system was operated in slower propeller driven types of 
aircraft, the system suffered erratic operation (intolerable errors), 
only when thunder storms set up exceedingly heavy corona conditions. Also, 
the system operated normally under conditions of moderate corona indicating 
that the low density types can be ignored. 

For the case of the jet aircraft, the corona discharge intensity would 
be mainly a function of speed and since the latter is usually close to 
maximum designed speed, the intensity will always be close to maximum, 
with corona discharge pulses occurring at all of the places where plastic 
surfaces are exposed to the airstream. The mathematical theory to des- 
cribe this phenomena is covered below and is closely allied to the more 
familiar treatment of fluctuation noise. 
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3* A development of Mathematical Description of Random Noise. 

Consider a plane-electrode diode for the emission limited case. A 
graph showing the current versus time for a single electron leaving the 
point of zero potential ( t - 0 ), and arriving at the point of highest 
potential ( t • t n ), is shown in Figure 3. When a constant potential is 
assumed over a time interval 2T, and a pulse having the shape of that 
shown in Figure 3> takes place for each electron crossing, a plot of 
the current l(t) versus time for the interval 2T is shown in Figure U. 

Assume that N electrons are emitted during the interval indicated. 
Then I(t) is the sum of the individual current pulses 




1-1 



and the time average is 




1-2 



substituting 





1-3 
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FIGURE 3 
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n 



Graph of i(t) versus Time 




FIGURE 4 

Graph of i(t) versus Time Showing the Interval 2T 
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Now consider the electron stream to be made up of a great number of 

19 

electrons whose fluctuations are considered to be random. A count of 
electrons arriving at the plate in the interval 2T is N. Repeating this 
count for a total number of M times gives a value for the average number 
of electrons per second as 



V = 



M 



- GO 



H, t- Nz->- Ms + • • • * ^ 

hA (at) 



i-U 



he M is increased and T held constant, some of the N's will have 
the same value, thus the number of N's of a particular value will in- 
crease as M is further increased. The probability of getting N electrons 
for some particular trial is defined as 

p(N) = lim Number of trials giving N electrons 

M-«-oo H J 

Since p(N) also depends on T it can be assumed that the randomness of 
the electron stream is such that the probability that an electron will 
arrive at the anode in the interval (t, t + i t) is 
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V where At is such that vAU<l, and that this probability is in- 
dependent of what happened before time t or what will happen after time 
t+At. For a very large number of trials this is a description of the 
Poisson distribution and is expressed as 

p(N) = {znf e- 2vT x - 6 

NJ 

Thus a large number, M, of intervals of length 2T having exactly 
N arrivals will be to a first approximation Mp(h). For a fixed value of 
t, equation 1-3 may be written 

%(t) r- 7 x - 7 

*J= » 4*2- > 

and for it effectively becomes 

Ut) - dt x _g 

\J~O0 

It is now possible to write the average l(t) over all of the M in- 
tervals instead of only over those having N arrivals, so as M 

co 

W = -1L p(N) %(t) 1-9 

N-o 
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Campbell's theorem 1 ^ states that the mean square value of the fluctuation 
about the average is 



(I(t) - I(t)) 



v V i 2 (t) dt 

>-o*> 



1-12 



i 2 it) - 2 i(t) i(t] - i(tr i-i3 



Since 



■>> 



■it 



i(t - tn) i(t - t m ) 



Y\i\ rr\s\ 

Averaged over all values of t^, t2> t^, 
fixed 



l-lli 



tj^ , with t held 



I>’IL X - 1S 



The multiple integral has two different values, so tfien n ■ n 




and for n / m its value is 




dL*»> 

■ZT 
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Counting up the number of terms in the double sum shows that there are 
N of them having the first value and N 2 - N having the second value. 



lii 



For 



- T-rk < t 



T - A 







i fy) dt ■/- 




Averaging over all the intervals instead of only those having N 
arrivals, the final form can be written as 



where the summation with respect to N is performed as in (l-ll) and after 
the summation, the value (1-11) for I(t) is used. 

The relationship of the above development to the problem can be shown 
for a small increment of surface area located at one of the -points where 



By equation 1-1 the individual current pulses can be described for a 
single element, and when all of the points are taken into consideration, 
plus all of the locations where plastic housings are adjacent to metallic 
airframe, equation 1-20 expresses the total average squared current. 
(Appendix I) 




1-19 




1-20 



the metallic airframe and plastic housing are adjacent. 
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1*. Previous Statistical Noise Measurements 

Some investigation of other statistical approaches to the measure- 
ment of noise was made before any design considerations were conceived. 

In their article, "A Statistical Approach to the Measurement of At- 
mospheric Noise" by R. S. Hoff and R. C. Johnson,^ it states in part,"... 
that atmospheric noise, causing noise voltage to be induced into a receiv- 
ing antenna. . .the envelope of this noise voltage after being amplified 
by the receiver stages and further modified by limiting and other non- 
linear processes, is an irregular function of time or time series l(t), 

111 

and may be treated by statistical processes." Also, Landon has shown 

that the distribution of noise voltage when passed through a narrow band 

filter is normal (Gaussian), (Appendix I) that is, the fluctuation noise 

output of a band pass amplifier has a distribution of voltage versus time 

which follows the normal error law. The method used is similar to that used 

l8 

by Lord Rayleigh in determining the probable distribution with time of 

the instantaneous amplitude of a large number of sine waves, all with 

12 

equal amplitudes and random relative phases. From Knudzton's report , "... 
that the narrow band noise has the character of an amplitude and frequency 
modulated wave for the ideal rectangular band pass filter, where the 
average period of the envelope is a linear function of the bandwidth only." 

Thus it is possible to list the conditions that will effect the de- 
sign characteristics of the device to examine the properties of the pre- 
cipitation noise. 

(1) If a quantitative value of the noise power is desired, a narrow 
band pass filter of approximately 3 kc must be used since the center 
frequency is 126 kc . (See Appendix II) 

(2) The device must simulate a navigational system which has the fol- 
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lowing characteristics: 

(a) Bandwidth 



30 kc 



(b) Resonant Frequency 126 kc 

(c) Servo time constant response - 100 seconds 

(3) Ideally considered, the noise has a white spectrum. This of course 
is not true, but is permissible for this problem since Plank's radiation 
formula indicates that the thermal noise energy is down l/e of its maxi- 
mum value of KT at about 10? me $ 

To satisfactorily solve this dilemma of having two bandwidths it would 
be necessary to have two distinct systems; the first to measure a quanti- 
tative value of the noise, and the second to ride next to the navigational 
system recording the relative noise input level for a specified naviga- 
tional track. It was decided to use the latter system to study the effects 
peak disturbance periods as a function of navigational error and to use 
the first system to set up a method to investigate the power spectrum of 
high density precipitation static by measuring the average power level 
across the band in small increments of frequency. 
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Chapter II 

Description of Designed Measuring Devices 

1. SYSTEM I - MEAN SQUARE NOISE METER 

As shown in Appendix II, it is possible to represent the output of 
a narrow band filter by 

e(t) = e cos (t) COS W m t + e sin ^ sin V 

Using a very narrow band pass amplifier as the input stage, the mean 
value of the output noise voltage can be found by integrating the voltage 
over a long period of time. This integrator has an output which is made 
up of a component that increases linearly with time, plus some variation 
due to the input a-c components wnich are smoothed out after sufficient 
elapsed time. It is specifically desired to record the mean square value 
of the noise voltage ( e^), by constructing a device which will square the 
input and then integrate the low frequency (difference terms) components 
over a long period of time. The advantages are: 

(a) The recorded measurements are proportional to energy. 

(b) From the basic premise, that the noise has a normal distribution, 

it must then consist of a uniform spectrum of Spectral Density N o 
2 

volts /cps at all frequencies other than zero, and have a true mean value 
(e)volts, it is necessary to continue the integration process for a length 
of time 

T = (2 x 10 h N q ) /P 2 I 2 sec 

in order to insure a 95% certainty that the measured value lies within 
P percent of the true mean.^ 
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Frequency Analysis of tne Square law Circuit: The current of a 



diode may be represented by the infinite series 

n 



• - 2 3 

!p “ c^e +* c 2 e -v- c^e - 



' ^ c n e 



in which the constants c n may be determined as are the coefficients in a 
Taylor series. 

Since the diodes are operated in a region where the dynamic character- 
istic is not a straight line, let it be assumed that the first two terms 
accurately represent the portion of the characteristic used. Thus 

*p = C 1 e + c 2 e2 



Substituting equation III-7: Appendix II 



1 p = C 1 ( cos "m 1 + e s (t) sin w m t) 

+ c 2 ( e c ^ cos w m t+e s ( t ^ sin w m t ) 2 

The only terms of interest are those containing difference terms since 
the integrator is essentially a low pass filter, therefore in this case 
it is 



2 c 2 e c (t) e g (t) sin w m t cos w m t 
By trigonometric identities this is equal to 

2 c 2 e c (t) e g (t) j^sin ( w m - w 0 ) t + cos w m t sin w 0 t^ 

Thus the presence of the difference terms allows the use of a large 
time constant in the integrator circuit to average the input low fre- 
quency terms. 
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Description of Circuits in System I 

Figure 1-0 Block Diagram of the System 

Figure 1-1 Narrow Band Pass Amplifier - Essentially a tuned grid, 
tuned plate tank, measured to be 3kc wide at half power points, with center 
frequency at 126 kc, and tunable over a 30 kc band. 

Figure 1-2 - Logarithmic Amplifier - The amplifier is an instan- 
taneous linear amplifier designed to obtain a large linear range. A feed- 
back loop from the Miller integrator circuit which has a time constant of 
100 seconds causes a compressing effect to the input signal. This allows 
examination of signals whose magnitude may vary from 10 micro-volts to 
10^ micro-volts. The variable gain is calibrated to give approximately 
10 db gain for - 6 volts of avc and 70 db gain for - 1 volt avc. 

Figure 1-3 Square Law Circuit - This circuit was originally designed 
using the relation 

E o = k e/ 

where E 0 is the output voltage and e| is the input voltage of the circuit. 
Knowing the desired avc characteristic and the time constant of the Miller 
integrator it is possible to approximate the value of k. This was accom- 
plished by plotting several curves using values close to the approximated 
value of k, constructing a circuit with variable resistors and adjustable 
bias levels, empirically determining the component values after the input 
and output values ' fit * the desired curve. 

Figure 1-U Detector and Miller Integrator - The detector is a peak 
reading device with a time constant determined by the center frequency 

T C detector = 1/1.26 x 10^ = 8 micro-seconds 
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The basic principle of the Miller integrator in this system is the 
operation of averaging. It is performed by passing the general input through 
a low pass filter (integrator), thus the contributions to the output will 
be made only by those terns which have frequencies which fall within the 
range of the filter. In this problem the observation interval is fixed 
at 100 seconds, so thus the output, M(T), of the averaging filter may be 
related to its input by means of the convolution integral involving the 
filter's impulse response and weighting function, h(t). 



bandwidth, and the weighting function, h(t), are given for an integrator 
as 



i 

M (T) = i(t) E (T - t) dt 




where T is the duration of the observation interval. 



From "Statistical Errors in Measurements''^, Table Ij the effective 



h(t) = l/T = o.Ol where 0 > t>T 

\ = 1/2T = 0.005 



The input will be the difference terra 



E d (t) 2 c 2 e c ( t ) e s (t) sin (w m - w Q )t 



for 



0 < w < 0.01 



Simplifying M(T) 
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The above represents the output of the averaging filter for an 
observation interval of 100 seconds. 

Figure 1-5 Meter Coupling Circuit - This circuit uses an Esterline 
Angus type recorder, Model AW. The meter movement is a 1 ma type. 
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SQUARE LAW CIRCUIT 
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2. System II - NAVIGATIONAL DEVICE SIMULATOR & RECORDER 

As was previously stated, this system must duplicate a navigational 
device that is represented by the diagram in Figure 1. It must record 
noise voltage but must disregard the signal, so that if the circuit shown 
in Figure 1 were used in its entirety, a recording of the output would 
simply be the average value of the input signal and noise lumped together 
with no means of distinguishing one from the other. It is required then, 
that the signal must be eliminated and only the noise voltage recorded. 
Filtering of course will be useless since it removes those noise frequencies 
that are of interest, so that the only solution is to phase out the signal. 
This will require two separate channels and the use of phase shifters, but 
in keeping with the basic premise, to duplicate the navigational system 
as much as possible, the only change of significance will be the intro- 
duction of a balanced modulator. Its purpose is to ict as a gate, opening 
and closing at a set rate, but 90° out of phase with the gate in the other 
channel. Thus when the two channels are added together, the signal will 
have been phased out but the noise components will be the sum of all the 
cosine terms (from one channel) and all of the sine terms (from the other 
channel) . 
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Description of Circuits in System II 
Figure II-O Block Diagram of System 

Figure II-l Band Pass Filter and Tuned Amplifier - Essentially tnis 
circuit is a single tuned, link coupled circuit, acting as a band pass 
amplifier with bandwidth of 30 kc with a resonant frequency of 126 kc . 
Figure H-2 Logarithmic Amplifier - Same as System I. 

Figure II-3 126 kc Crystal Oscillator and Phasing Networks. A highly 
stable Pierce crystal controlled oscillator is used as a local oscillator, 
the output of which is fed to two RC type phasing networks. If the out- 
put of the oscillator is 

e Q ■ E sin w Q t 

then the outputs of the phasing networks will be respectively 

■ E sin w Q t - U5° 



and 



e 2 = E sin w Q t ■+■ 1£° 

Sin&e e^ and e ^ are ninety degrees apart and bear no fixed relationship to 

the input signal as yet, they can be represented as 

e-i = E sin w t 
-L 0 

e 2 = E cos w Q t 

Each are amplified in separate channels and remain in these channels un- 
til detections take place. 

Figure II-h 126 kc Amplifier, Gate and Integrator - The first two stages 
are voltage gain stages which split e 2 into two voltages 180° apart in 
phase and thus drive a push pull cathode follower stage. By this means 
it is possible to offer a low impedance driving source to the balanced 
mo dulat or ( ga t e ) . 



30 



Frequency Analysis of Coherent Detector (Balanced Modulator) 

If the input noise voltage has the form at the output of the log 
amplifier of 

e(t) a (t) sin wt *♦" A2 (t) cos wt 

where 

w = 2 It ( 126 kc * 30 kc ) 

the output of the balanced modulator can be computed in the following 
manner. 

Assuming matched diodes with similar current characteristics, the 
current can be represented by the infinite series 

l 

S> = a o 4 a i e a 2 e 2 +■ a 3 e 3 ....+ a n e 11 
Since the integrator will act as a low pass filter on the difference 
terms of the output of the non-linear device will be of consequence. 

The d-c term (a Q ) can be ignored since transformer coupling is used. 

The voltage at point (1), Figure 11-5 is 
Ei = e^e(t) 

The voltage at point (2), Figure II-5 is 
E 2 = - e(t) 

thus the value of e can be written as 

out 

A(E 2 ' e o) 2 “ A(e 0 - F^) 2 
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Setting AR - 1/K 



K e o = (" e i' e e U) “ e 0 ) 2 " (e D - e x - e(t)) 2 

Thus 

- U e, e(t) 

e o s I — 

K - ii e x 

Substituting 

• It E Ai (t) sin w t sin wt - U E A 0 (t) sin w t cos wt 

6 2 O c. Q 

K - U E sin w t 
o 

Since only the difference terms are of consequence 

- sin W Q t cos wt * - sin (w Q - w) t - cos w Q t sin wt 

and 

sin w Q t sin wt - - cos (w q - w) t cos w Q cos wt 
Thus e Q can be represented as a voltage having frequencies made up of 
the difference terras. 

For channel two the analysis is the same with the exception of the 
value of eg 

-U e 9 e(t) 
e o s 2 . - 

K - U e 2 
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Substituting 

- -1* E Aj^ (t) cos w Q t sin wt - l* E A g (t) coS w 0 t cos wt 

K - I* E cos w t 
o 

Again only the difference terms are of consequence 

- cos w Q t sin wt * sin (w Q - v)t - sin w^t cos wt 

and 

- cos w Q t cos wt * sin w Q t sin wt - cos (w Q - w) t 

as before, the voltage applied to the integrator can be represented as 
having frequencies made up of the difference terms. 

The Miller integrator is similar to that used in System I except 
that the output can be either a positive or negative value, since the 
average will vary about zero and will have a mean average level of zero. 
Figure II-5 1*00 Cycle Oscillator and Amplifiers - The 1*00 cycle os- 
cillator is a Colpitts type with output fed to the two channels after 
being amplified and shifted in phase 90° with respect to each other by 
the two RC type networks. The method is similar to that previously men- 
tioned in connection with the 126 kc oscillator and the shift is in keeping 
with the previous phase shifts. The purpose is to obtain the sum of the 
two channels across a transformer. It should be mentioned that the- shift 
to 1*00 cycles is to enable a power gain to be obtained before taking the 
sum of the two channels. 

Figure II -6 1*00 Cycle Chopper and Amplification Stages. - The 1*00 cycle 

chopper is a Model manufactured by the Bristol Co. of Waterbury, Conn.. 
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The necessary power amplification is provided by a single stape in each 
channel with the two channels being coupled together using a push-pull 
amplifier. The output to the detector is the sum of the two channels. 
Figure II-7 Detector and Recording Circuit - The detector is a peak 
reading type with a time constant proportional to the I 4 OO cycle fre- 
quency. The Recorder is an Ester line Angus, Model AW, having a 1 ma 
movement. 
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LOG AMPLIFIER 
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126 KC AMPLIFIER, GATE a MILLER INTEGRATOR 
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DETECTOR 8 RECORDER CIRCUIT 



3. OPERATION 



The technique of taking data is quite simple. 

Previous calibration of the meter and circuits should be made in 
the laboratory using a General Radio Signal Generator 803C , fixed at 126 kc 
with zero modulation. To simulate an antenna load a value of capacitance 
equal to the antenna capacitance of the particular aircraft model should 
be shunted across the input terminals so that the tuning capacitor can be 
set for the proper value. 

For an input of ten microvolts the meter should read a zero value 
while for a full scale reading a value of 100,000 microvolts can be used. 
Intermediate points can be located to suit the user. 

Flights must be at least four hours long in order to accumulate suf- 
ficient data per flight. It has been found that a flight of this duration 
is necessary mainly to show the magnitude of the error of the navigational 
device for a known course and speed over a closed geographical loop of known 
dimensions . 

Flights in the vicinity of thunder clouds are desired with a record of 
the time spent in such locality marked on the recorder paper along side of 
the measurement record. 



CHAPTER III 



RESULTS AND CONCLUSIONS 



1. Results: 

System II was calibrated in the laboratory and will be delivered to 
the Air Force some time in April 1 9$k- No particular recommendations or 
conclusions can be offered at this time since no flights have as yet been 
made, but as is typical with statistical data, the majority of time spent 
in investigation is used to gather data. 

System I was constructed and calibrated in the laboratory during the 
month of April 19$h» A unique testing set up, developed by Tanner^ was used 
to generate corona pulses and is illustrated in Figure 5>. 

Preliminary recordings of the output of this test setup indicated only 
minute amounts of power were being generated. It was decided to take the 
unit on a test flight to ascertain whether any general indications of pre- 
cipitation static could be obtained. 

A flight was made on April lit, 19^U from 0830 to 1230 PST with route 
as follows: 

From Monterey, California to Fresno, California, and then north-east 
to the Sierra Mountain chain, following this due north to Yosemite National 
Park. 

A long line of cumulus type clouds hung over the ridges offering good 
precipitation static conditions, so the system was placed in operation and 
immediately the recorder indicated moderate corona because all discharge 
wicks were removed from the plane previous to take-off. Figures III and 
IV are parts of the recordings taken during this flight and show vividly 















. 






SYSTEM I 
RECEIVER 




FIGURE 5 



Corona Pulse Generator 
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the difference in the average power levels for an area of negligible pre- 
cipitation static and an area of heavy precipitation static. 

Conclusions: 

Investigations of precipitation static made using System II will 
record values related only to the average power input to a particular 
device so it is obvious that this information is of no value to anyone 
interested in the actual power spectrum or the possible frequency distri- 
butions of precipitation static. It was for this purpose that System I 
was designed. 

Since the data gathering period will take considerable time, only 
a brief outline of the method and approach will be given. 

The approach to the problem can be again allied to the familiar treat- 
ment of fluctuation noise, using the autocorrelation function and its 
relation to the Fourier cosine integral. 

The autocorrelation function / (Y) of a time function i(t) is defined 
as the time average of the product of all pairs of points i(t) separated 
by an interval n . a T~ , 

96 W tfc) L (t +T) 

Substituting from Equation 1-20 — 

/(nr) - v c (t+ r) <Lt t T(i) x 

In the above the second term is the square of the d-c value I (t). 

If a current i(t) is defined as the fluctuation of the current about its 
mean value, then 
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i (t) = i(t) - I 

pOO 

fAr) = v \ i(t) i (t + t ) 

± J-0O 



Note that the total’ noise power depends not only upon the average current 
but on the shape of the current pulse as well. 

The noise power spectrum G^(w) can be obtained from the autocorre- 
lation function by the Fourier cosine integral relation, and is found to 
be 



where e is equal to the charge of electron, and can be plotted as shown 
in Figure 6. 

The noise power in the frequency range f^ to fg is 



In the frequency range for which wT*.^*^ , G(f) is nearly constant 
in value. This is true no matter what the shape of the current pulse is 
since the value of the Fourier Integral at zero frequency depends only 
upon the zero moment of the function. Thus G(0) depends only upon the 
area under (r) and not upon its shape, so that the noise power in a 
bandwidth A f is 2 v eA fo 



G i (w) = 2 v 



^ 6 (w ToS)^ 2(1 - cos wT*,- w sin w T^) 

( W ToJ U 



^G(f) df 
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FIGURE 6 

Plot of Noise Power Spectrum versus oj*! 
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A device then, having a bandwidth that is less than one radian can di- 
rectly measure increments of the power spectrum. This work is now in pro- 
gress using System I using the noise model described in the beginning of 
this chapter. It is obvious that the noise model device is inadequate as 
a generator for heavy corona conditions and investigation and design of a 
superior device is contemplated. 
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APPENDIX I 



THE NORMAL DISTRIBUTION 

1. Probability: 

(a) If an event is certain to occur it is assigned an ar- 
bitrary value of unity. For example: 

In the toss of a single die the probability of 



a six 


is 


1/6 


a five 


ii 


1/6 


a four 


1! 


1/6 


a three 


11 


1/6 


a two 


11 


1/6 


a one 


11 


1/6 



thus, if six probabilities are mutually exclusive, the probability that 
either one or the others will occur is equal to the sum of the separate 
probabilities. 

(b) When a particular result in any trial is independent of 
the results obtained in any other trials, is called a random process. For 
instance: 

A large number of independent trials in tossing a coin mi^it reveal a 
consecutive string of heads. This will not affect the next toss since the 
probability is still 1/2 that the coin will come up either a head or a tail. 

(c) The probability of the simultaneous occurrence of n events 
is the product of the n individual probabilities of n events. Thus the 
probability of winning all tosses in a season of eight footballs games 

O 

would be (0.5) = 0.0039 

The above three rules form the basis of the probability distribution 
laws upon which the mathematical description of noise is based. 
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2. Distribution formulas: 



If the symbol P m (n) is made to stand far the probability that an 
event will occur exactly n times in m independent trials, then the basic 
distribution formula, (Bernoulli Distribution), is 



P (n) = 

m' ' 



ro n, v m-n 
c n P (1 - p) 



n m-n 

— tV r-r P (1 " p) 

n I (m - nj I 



where p is the probability for the event to occur in any particular trial. 
This specializes to, (Poisson Distribution), 



for m » n 

p ^ 1 

m ■?> 1 



P(n) = -n 
7 



(n) n 

~ ill 



where 



n 



- mp 



When n becomes very large as in the case of hi^i density Poisson 
noise, the above further specializes to the Gaussian Distribution. 



P 

m 



(n) 



when 



1 e -(n-n) 2 / 2n 

-fT#TT 

n >> 1 



n-n 



n 






1 



n - H > > l 
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3. The Central Limit Theorem 

The Central Limit Theorem states that subject to very lenient 
conditions, and regardless of the amplitude distributions of the indi- 
vidual components, as the number of components n approaches infinity, 
the amplitude distribution of the current approaches the Gaussian or 
normal distribution. 

Thus it is possible to relate the limiting form of the Gaussian 

Distribution to the Mathematical Description of Noise (Chapterl) by 

virtue of the large number of independent events (electron crossings) 

✓ 

which add up to give the independent variable I. This distribution 
will be independent of the amplitude distributions of the individual 
events (shapes of the electron pulses). It is required, generally speak- 
ing, that the contributions of the individual components be small compared 
with the total value of the sum 



APPENDIX II 



PROPERTIES OF NARROW BAND NOISE 

1. Examining the statistical properties of the output of a narrow 
band filter plotted in Figure A-I, and the assumption, that the band 
width of the filter, f m ^ jthe output can be expressed in the form 

e(t) s E(t) cos (w m t+^(t)) II-l 

provided that the input noise has a Gaussian distribution. 

E(t) then represents the magnitude of the envelope of the output and 
^(t) the phases changes. 

Using the Fourier representation of the output as 

e(t)= / (a cos nw t + b sin nw t) II-2 

n o n o 



where the summation is to include all frequencies which is in the range 
of the filter. 

Setting up the sum and difference terms so as to show relationship 



to a sine wave 


of 


frequency w m j' 






COS 


nw t 
0 


= 


cos ^ 


L (n "° 


-w )t +■ w t ' 
n m _ 


l 


II-3 






= 


cos 


(nw 

o 


-w m )t cos w m t 


- sin (nw -w 

o m 


) sin w t 
ra 


sin 


nw c t 


= 


sin | 


L (nW o 


-w )t 4- V t 1 
m raj 




ii-U 
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sin 


(nw Q 


- w)t cos w t + cos (nw 
m m ' o 


- w m )t sin w B t 
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FIGURE A- I 

Output of a Narrow Band Filter 
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Letting e cos (t) represent the summation of all the coefficients of 

e°s w m t in the expression for e(t), and ©gi^t) represent the summation 

of the coefficients of sin w t. then 

m ’ 



e cos (t) s I [fn cos (nw o - V t + b „ (<“ 0 - w m )t] 1 1-5 

and 

A - — p - 

e sin (t) = V_ b n cos “ w m )t - a n sin (nw Q - w m )t J II-6 

v\ 

Substituting 

e(t) = e cos (t) 008 V + e si„ V n ‘ 7 

Thus, e (t) and e . (t), are slowly varying compared with 

COS 81.11 

cos w t or sin w„,t. 
m m 

From the Expression III-l, the magnitude 



and 



E(t) = 



1 



COS 



(t) +■ e 



sin 



(t) 



II-8 



(t) - tan 



-e . (t) / e (t) 

sin cos v ' 



II-9 



The random variables e (t) and e . (t) are sums of independent 

cos sin 

variables with Gaussian amplitude distributions. The average power can 
then be written as 
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E 

2 



11-10 



cos 



2 

e . 
sin 



- e 
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'cos 



»I -5 



n 
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sin n n 
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